% clc;
% syms a b positive,
% syms x;
% f=a*sin(8*x^2)+b*cos(-2*x+2);
% L=feval(symengine,'limit',f,'x=infinity','Intervals')

syms a b positive, syms x; f=a*sin(8*x^2)+b*cos(-2*x+2); L=feval(symengine,'limit',f,'x=infinity','Intervals')